Interaction-assisted SU(2) adiabatic holonomies in Josephson phase qubits 
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We propose a scheme for generating SU(2) adiabatic geometric phases in a circuit consisting of 
three capacitively coupled flux-biased Josephson phase qubits. 

PACS numbers: 03.67.Lx, 03.65.Vf, 85.25.-j 



a^ ■ 
o : 
o . 

(N ■ 

5-H ! 
Oh. 

<\ 

^ : 



^ ■ 

^ ; 



> 

00 
(N 
O 

o 

O 



X 



I. INTRODUCTION 

When slowly changing, external controls that gov- 
ern adiabatic evolution of a quantum system return to 
their original values, the phase of system's wave func- 
tion acquires a purely geometric contribution indepen- 
dent of both system's energy and duration of the adia- 
batic process. This U{1) contribution, initially discov- 
ered by Berry [ij, mathematically interpreted by Simon 
0, and later generalized to non-abelian U{N) operations 
by Wilczek and Zee |^] , has been the source of many in- 
teresting developments in various fields of chemistry and 
physics [3], such as molecular and spin dynamics, op- 
tics, fractional quantum Hall effect, various branches of 
condensed matter physics [5| , and others, such as, for ex- 
ample, the field of quantum information processing 

Of particular interest to us is the possibility of observ- 
ing non-abelian geometric phases (SU(2), in our case) 
in superconducting circuits with Josephson junctions 
Such circuits are currently considered as promising can- 
didates for scalable quantum computing architectures [1] . 
Even though in this paper we consider a definite phys- 
ical system (the so-called capacitively coupled flux bi- 
ased Josephson phase qubits 9]), our analysis may ap- 
ply equally well to any quantum computing architecture 
whose dynamics is described by the Hamiltonian given 
in Eq. ([7]). This very simple Hamiltonian has properties 
that make the corresponding system ideally suitable for 
actual realization of non-abelian geometric phases. The 
most important of these is the presence of two doubly 
degenerate subspaces, one of which, Vi, can support non- 
abelian phase, and the other, V2, cannot. Experimentally 
then, V2 may be used as fiducial, reference subspace rela- 
tive to which various interferometric experiments involv- 
ing Vi could be performed. 

To recall how geometric phases appear in quan- 
tum mechanics and to fix our notation, let us con- 
sider a quantum system whose Hamiltonian depends on 
the time t via several controllable parameters X{t) = 
{X^{t), X'^it), X^{t), ...},n^ T-l-iHt))- It is assumed that 
H{X{t)) forms an iso-degenerate family of Hamiltonians 
(no level crossings as A varies) @ . Let the instantaneous 
energy eigenbasis Xa{X{t)) within a give degenerate sub- 



space be chosen, 

n{Xit))xa{X{t)) - E{X{t))xa{X{t)), a = 1, 2, . . . , TV. 

. 

At t = 0, prepare the system in one of the energy eigen- 
statcs, ipa{0) = Xa{oi{0)), and let it evolve according to 
the Schrodinger equation, 



ld2ljait)/dt^n{X{t))ija{t). 



(2) 



Then, in accordance with the adiabatic theorem, at some 
later time t, the state ipa{t) will be a linear superposition 
of various Xb{X{t)) belonging to the same degenerate sub- 
space. Therefore, in general, 

Mt)^e-^foEiHt'))dt'Y^^^^X{t))UUm), (3) 

b 

where the first factor on the right is the usual dynamical 
phase, and U{X{t)) is the unitary matrix representing 
additional, purely geometric "rotation." The U{X{t)) is 
the sought for non-abelian geometric phase. It is given 
by the path-ordered exponential, 



U{X{t)) =7'exp. 



^*^A(A(0) AXO^ij, (4) 



where 



[A{X)]cb^{xc{\)\d/dy\xb{X)) 



(5) 



is the so-called adiabatic connection, also known as the 
M(iV)-valued gauge potential. In general, when X{t) 
traces a closed loop C in the parameter space, the system 
picks up a nontrivial geometric phase (the holonomy) , 



Uc 



7'exp|-^^A(A)dA*| ^ 0, 



(6) 



independent of the speed with which the loop had actu- 
ally been traversed (as long as it was done adiabatically) . 
In the next section we describe a multi-qubit solid state 
quantum computing architecture capable of supporting 
such a phase. 



II. CAPACITIVELY COUPLED TRIPARTITE 
SYSTEM AND ITS HOLONOMIES 
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One system that leads to nontrivial, SU(2) holonomies 
consists of three capacitively coupled Josephson phase 



2 



qubits [Ol whose effective Hamiltonian in the rotating 
wave approximation is given by (cf. [10]) 

n{B.,,By) = (l/2)[B,K+4)+B,K + (Tf) 

+9 Kf^2 + ^3^2 + ^3^f + ^l^l) 

+J{a-af+ay<jy)], (7) 

where, in polar coordinates, Br^ — B cos (f>, By — B sin (p. 
Here, qubits 1 and 2 experience the same magnetic field 
produced by identical microwave driving on both qubits. 
This microwave drive is the parameter A = {B, (j)} that 
will undergo adiabatic change. For future use we define 

A = (1/2)(-B2 + 6g')/^B^ + Ag\ 

A' = (1/2)(B2 + 6g^)/^Bi + Ag^ > 0, 

K = B^ - 2g'^B^ + 8/ > 0, 

L={-B'^+ V)^B4 + 454. (8) 

The corresponding Hamiltonian matrix 'H{B,(f)) has 
six different eigenvalues, two of which are doubly degen- 
erate. The non-degenerate ones are Ei^fij^s — (l/2)[>/± 
■\/ + 4{2g^ + B'^ zL 2Bg)]. Their corresponding eigen- 
kets are rather complicated, each leading to the same 
C/(l) Berry phase e*'''^ = e"'^'^* on a single precession of 
the field, and will not be needed in what follows. The de- 
generate eigenvalues that will be needed are £"3^4 — — J, 
with 

and £^1^2 = 0, with 

= {{-B^+4g^)\000) + B^e^'^\0n) 
-2gBe"f'\l00)}/V2K, 
|X2(S,<^)) = {25i?3e-3^'^|000)-4g3Be-^'A|011) 
-B^{B^ ~2g^)e-^"l'\l00) 

+K\lll)}/^/2K{BU~ig^. (10) 

A straightforward calculation based on Eq. ^ then gives 
the adiabatic connection within the |xi,2) subspace, 

dA{B,(j)) = i{[2B\B^ + 2g^)/K]{a'-/2) 

-{2gB^A/K) (cr^ cos 30 -I- cr^ sin3(/>)} d<j) 
-i{2gB^A'/K){a'= sinSt/)- ci^ cos3(/)) dB. 

(11) 

Here, the Pauli matrices operate within the lxi.2) "ge- 
ometric" subspace and are different from the sigma- 
matrices originally used to describe the qubits in the cir- 
cuit. 

Notice that in the course of adiabatic evolution we are 
not allowed to take any of the limits J — > 0, 5 ^ 0, 
or i? 2g, since each breaks the requirement of iso- 
degeneracy. Also notice that J appears nowhere in |xi,2), 
£'1,2, or dA, so the geometric phase is not sensitive to the 



precise value of J. Nevertheless, the presence of J 7^ is 
crucial here: it guarantees separation of the |xi,2) sub- 
space from the rest of the Hilbert space. Such (at least, 
partial) robustness against possible imperfections in the 
coupling is important since our ability to implement mul- 
tidimensional geometric phases relies on the exact de- 
generacy of the underlying energy subspace. However, 
in order to have the exact degeneracy of i?i,2, we also 
require that the coupling g be the same for 1 3 and 
2 3. This is a very important requirement, not eas- 
ily achievable experimentally (unless tunable coupling is 
available). Additionally, since B does appear explicitly 
in Eq. (|10[) . we require good control over B, which is not 
a problem experimentally. In a typical experiment, even 
for long times, the magnetic field (the microwave drive) 
can easily be controlled [11]. 

We also point out that the sequence of operations gen- 
erating non-abelian phases must be short compared to 
the decoherence and dephasing times. On the other hand, 
system's evolution must remain adiabatic compared to 
the couplings g and J and the magnetic field B. In 
a typical experiment, g is on the order of 20 MHz, or 
1 /g ~ 50 ns, so the entire operation would be several hun- 
dred nanoseconds which is comparable to the presently 
achievable decoherence times. However, with stronger 
couplings up to 5 = 100 MHz and 1/5 = 10 ns it is pos- 
sible to perform an adiabatic operation in several tens of 
ns [llj. 

Let us now fix a basis \xi,2iBo)) at = and initialize 
the system in some state iV'o) = o:\xi{Bo)) + P\X2{Bo)). 
We may then trace a closed loop C in the parameter space 
by performing the following sequence of operations: 

(1) Begin, if needed, by adiabatically changing the field 
from Bq to Bi along the radial direction in the (x, y)- 
planc (no precession). The system will pick up a phase, 

e{Bo,Bi) = 2 [ \2gB^A'/K)dB. (12) 

J Bo 

(2) Perform adiabatic precession at Bi, 

Bx{t) — Bi cos Wot, By{t) — Bi sin Wot. (13) 

If the field makes n revolutions of total duration = 
2Tm/Lj{), the accumulated phase will be 

Uxz = e_^^^{cos(2A7m) 

+isin(2A^n) [(25SVif)CT" + (L/if)a"] }|b=Bi , 

(14) 

which represents a rotation by angle a„ = 4ATrn\B=Bi 
about the axis n = {-l)"'-^[2gB^/K,0,L/K]\B=Bi that 
lies in the {x, z)-plane of the Bloch sphere (constructed 
with respect to |xi,2))- The set of gates that can be 
generated this way is shown in Fig. [T] 
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FIG. 1: The one-parameter family of gates in the \xi,2{B)) basis that can be generated in a single precession using non-abelian 
holonomy Uxz(B) of Eq. (I14p at n = 1. Panel (a): _B-dependence of the rotation angle. Panel (b): _B-dependence of the 
rotation axis. To implement a given gate (such as a Hadamard, NOT, PHASE, etc.), a specially designed computational basis, 
different from |xi,2(B)), must be chosen. See Sec. Illll for relevant examples. 



(3) Bring the field back to its original value. This wiU 
result in additional contribution to the phase, {Uy)^^ ■ 
The full holonomy is then 



Uc = {Uvr^UxzUy = UxziUyf 



(15) 



III. DESIGNING GEOMETRIC GATES 

To demonstrate how this non-abelian scheme works 
in simple situations, let us find a basis in which a use- 
ful gate, such as a Hadamard, can be made in a sin- 
gle revolution of the effective magnetic field. We want 
Uy = ^, n = 1, and Uq = Uxz = H, possibly up to 
an insignificant U(l) factor. First, notice that the choice 
I i> T) = lxi,2(-B)) does not work as can be seen from Fig. 
[TJ We therefore introduce a new basis. 



i> := cos/3|xi(B))+sin/?|x2(S)), 
T) := -sm(3\xi{B))+cosl3\x2{B)), 



(16) 



with some unknown, _B-dependent angle /3. Setting n = 1 
and changing to the new basis gives 

Uxz{B,(3) = - cos(2A7r) - i sin(2^7r) X 

{[{L/K) cos2(3+{2gB^/K)sm2P]a'' 
- [{L/K) sin 2/3 - (2gB^/K) cos 2/3]cr"^}. 

(17) 



It is then straightforward to check that by choosing 



B = 



/2{-12 + ^[39 - 4mi(mi + 1)](1 + 2mi)2} 



4mi(mi -|- 1) — 3 



f3H 



arctan 



m27r 



1 -I- {2gByL) + v/2[l + (2.gB3/L)2] 



1 - {2gByL) 



(18) 



with mi e {0,1,2}, and thus B/g « 1.9587, 1.3716, 
0.8375, and G ^, we get a Hadamard gate, Uxz = 
{zLi)H. On the other hand, by choosing same B and 
setting /3not = Ph — tt/S, we get another basis, in 
which Uxz = (±i)NOT. Other important gates, such 
as V i NOT and various PHASE gates, can also be gen- 
erated in a similar manner. 

We emphasize that the specially designed basis is only 
"special" in the sense that it allows a given gate (say, a 
Hadamard) to be generated by following a particularly 
simple path in the parameter space — here, by making 
a single turn of the adiabatically precessing in the {x, y)- 
plane field. 



IV. CONCLUSION 

In summary, we have shown how capacitively cou- 
pled tripartite superconducting qubit system may sup- 
port non-commuting SU(2) geometric holonomies. The 
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non-abclian character of the holonomies is due to the 
inter-qubit couphng. If we treat the relevant degener- 
ate subspace as a "geometric" qubit, then, with a proper 
choice of the computational basis, any of the standard 
single-qubit gates used in quantum information process- 
ing can be generated by the adiabatic transport alone. 



0404031. The author thanks Ken Brown, Michael Geller, 
and Matthew Neeley for helpful discussions. 



Acknowledgments 

This work was supported by lARPA under grant 
W911NF-08-1-0336 and by the NSF under grant CMS- 



[1] M. V. Berry, Proc. Roy. Soc. London A 392, 45 (1984). 
[2] B. Simon, Phys. Rev. Lett. 51, 2167 (1983). 
[3] F. Wilczek and A. Zee, Phys. Rev. Lett. 52, 2111 (1984). 
[4] See, e.g., A. Tomita, R. Y. Chiao, Phys. Rev. Lett. 57, 

937 (1986); C. A. Mead, Rev. Mod. Phys. 64, 51 (1992); 

D. Arovas, J. R. Schrioffer, and F. Wilczek, Phys. Rev. 

Lett. 53, 722 (1984); H. Koizumi, T. Hotta, Y. Takada, 

Phys. Rev. Lett. 80, 4518 (1998). 
[5] A. Shapere, F. Wilczek, Geometric Phases in Physics, 

World Scientific (1989); A. Bohm, A. Mostafazadeh, H. 

Koizumi, Q. Niu, J. Zwanziger, The Geometric Phase in 

Quantum Systems, Springer (2003). 
[6] P. Zanardi, M. Rasetti, Phys. Lett. A 264, 94 (1999); J. 

Pachos, P. Zanardi, Int. J. Mod. Phys. B 15, 1257 (2001); 

A. Eckert, et. al, J. Mod. Opt. 47, 2501 (2000). 



[7] L. Faoro, J. Siewert, and R. Fazio, Phys. Rev. Lett. 90, 
028301 (2003); M. Cholascinski, Phys. Rev. Lett. 94, 
067004 (2005); V. Brosco, R. Fazio, F. W. J. Hekking, 
and A. Joye, Phys. Rev. Lett. 100, 027002 (2008). 

[8] See, e.g., J. Q. You, F. Nori, Phys. Today 58 (11), 42 
(2005); G. Wendin and V. S. Shumeiko, Low Temp. Phys. 
33, 724 (2007). 

[9] M. Stcffcu, M. Ansrnann, R. C. Bialczak, N. Katz, E. 
Luccro, R. McDermott, M. Neeley, E. M. Weig, A. N. 
Cleland, J. M. Martinis, Science 313, 1423 (2006). 
[10] A. Galiautdinov, J. M. Martinis, Phys. Rev. A 78, 

010305(R) (2008). 
[11] M. Neeley (private communication). 



